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Abstract - In this paper we show how the stochastic quantization method developed by Parisi 
and Wu can be used to obtain Casimir forces. Both quantum and thermal fluctuations are taken 
into account by a Langevin equation for the held. The method allows the Casimir force to be 
obtained directly, derived from the stress tensor instead of the free energy. It only requires the 
spectral decomposition of the Laplacian operator in the given geometry. The formalism provides 
also an expression for the fluctuations of the force. As an application we compute the Casimir 
force on the plates of a flnite piston of arbitrary cross section. Fluctuations of the force are also 
directly obtained, and it is shown that, in the piston case, the variance of the force is twice the 
force squared. 



Fluctuation-induced (Casimir) forces [lip] are currently developed by Parisi and Wu 12 . These authors construct 



ttracting renewed attention, probably because of the ac- 
cessabiUty of small systems at the nano- and microscale, 
as it is in these small systems where Casimir forces are 
revealed, despite some macroscopic effects ^3j. However, 
(^techniques to calculate forces for complicated geometries, 
CO beyond the usual ones with high degrees of symmetry, are 
scarce. Recently, a powerful technique to calculate elec- 
tromagnetic (EM) Casimir forces [4] has been proposed. It 
is based on a multiscattering technique and has been suc- 
cessfuUy applied to many configurations, such as plates, 
. .cyhnders, spheres, wedges, etc. [s] 

Casimir forces have their origin in fluctuations of EM 
k>( fields. Already Lifshitz in 1956 ^ use the concept of fluc- 
5_J tuations (in this case of the electric currents inside a metal) 
C^to compute such forces. With the idea of fluctuations in 
mind, some authors [7|[8] used a method that takes as its 
starting point the Langevin equation that describes the 
evolution of systems subjected to fluctuations (in their 
case of thermal origin). Such method was applied to equi- 
librium situations, but also to non-equilibrium ones [9 



a Langevin equation for a given field subjected to thermal 
fluctuations, in a non-physical, pseudo-time. Despite the 
classical nature of fluctuations we will show how quantum 
fluctuations arise naturally within that method. We will 
follow the formalism developed in Ref. 11 , that takes as 



11 



a starting point the Langevin equation derived by Parisi 
and Wu, and construct its solution by diagonalization of 
the deterministic evolution operator with the appropriate 
bounday conditions — the immersed bodies — , and later 
averaging over the thermal noise. 

The resulting method has the advantage of computing 
the Casimir force directly and can be directly applied to 
compute torques and stresses over extended bodies. Also, 
it is suitable for numerical or perturbation methods, al- 
lowing to compute forces on realistic geometries. Finally, 
the effects of finite temperature are included directly. 

The Casimir force is calculated via the stress tensor, T 
(see Ref. llj), which must be averaged over the probabil- 
ity distribution of the fields 4>. Such a probability distri- 
bution is given by: 



in which the Fluctuation-Disspiation theorem is not sat- 
isfied. However, application of the Langevin approach to 
EM Casimir forces of quantum origin has not yet been 
developed. 

The goal of this paper is to apply the Langevin equa- 
tion method to calculate forces of quantum origin. The 
method is based on the stochastic quantization approach 
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-S[4>]/h 



(1) 



where S'[(/)] is the action of the scalar field with zero mass, 
Wick-rotated in the time variable (2| {t — ir), i.e.. 



S[(t>] 



dr / dr 



+ V^U, (2) 
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where /3 = l/kgT and Z is the partition function 
Z = J Dcjye-^'^'t'^/''. For the bosonic case, the field 6 
must obey periodic boundary conditions in time, that is, 
(j){T + f3h, r) = 0(t, r). The stress tensor is normahy a bi- 
hnear form in the field T = Ticf), 0, r] , the expression that 
defines the stress tensor operator T- 
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Herein we use the formalism of Parisi and Wu 
evaluate the average of the stress tensor in an alternative 
way. The idea of Parisi and Wu consists in a formulation 
of quantum mechanics or quantum field theory in terms 
of a stochastic process. More precisely, a Langevin equa- 
tion for the field (f) is written as an evolution equation 
in an auxiliary time, which we will call pseudo-time, s. 
In this description, the field depends on the new pseudo- 
time variable: (/)(t, r) (j){T, r; s). The Langevin equation 
takes the form 



is to reproduce the full quantum-thermal statistical prop- 
erties of the field. In particular, the complete quantum 
correlations, including quantum coherence, are obtained. 
The second difference is that the Parisi- Wu approach does 
not consider the fluctuations in the macroscopic metal but 
on the (electromagnetic) field itself. 

Having an expression for the stress tensor and a 
Langevin equation for the field, we can follow the proce- 
dure recently developed in to obtain the Casimir force. 
The field is written as (and a similar decomposition for 
the noise ij with coefficients ?7„m): 



a^(^(T,r;s) 



5(t> 



+ r]iT,r;s) 

+ wA(b + T^{T,v;s). (3) 



The term 7y(r, r; s) is the source of fluctuations, given 
by a Gaussian white noise satisfying the fluctuation- 
dissipation relation [m] 

{V{r,r;s)) = 0, (4) 
{T]{T,r;s)r]{T',r';s')) = 2kBTSiT - T')S{r - r')S{s - s'). 

Higher correlation of odd number of ?7-functions vanish, 
while correlations of an even number of ?7-functions fac- 



torize in all possible products of two ry-functions 15 . For 
instance, the 4-correlation: 

(77(l)7?(2),y(3)ry(4)) - i2kBTf [6(1 - 2)6(3 - 4) 

+S{l-3)S{2-4) + S{l-A)5{2-3)], (5) 

where the symbol 1, 2, 3 or 4 stands for the complete set 
r, r, s, and the S means the product of 3 (5-functions, as in 
Eq. ^. 

The solution of the Langevin equation in the stationary 
limit s — > oo reproduces the probability distribution given 
by Eq. ([!]) [l3] . The pseudo-time s has not physical mean- 
ing and the role of the evolution in s in solely to sample 
over the probability distribution ([I]). Statistical correla- 
tions of two or more fields, must be evaluated at the same 
value of s, but otherwise the positions and times r can be 
arbitrary. 

This approach, although similar, is qualitatively differ- 
ent to the Lifshitz method. In Lifshitz method, fluctu- 
ating electric currents are introduced in the conductors 
which produce the fluctuating electromagnetic field in the 
vacuum. The fluctuations of the current are thermal, tak- 
ing place in real time, and the quantum character of the 
model is introduced by the use of the quantum fluctuation- 
dissipation relation . In the Parisi- Wu method, the fluc- 
tuation process takes place in a pseudo-time and its role 
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0nm(s).9m(T)/„(r), 



(6) 



where /„ and gm are the eigenfunctions: 
V2/„(r) = -A2/„(r), l|^g„(r) = 



"W™5m(r), (7) 



which are orthogonal under the L^ scalar product in space 
or time. The above expression indicates that /n(r) and 
encode the spatial configuration of the system, that is, 
the position of the bodies and their boundary conditions. 
In a similar fashion, gm{T) and uS^ contain the (Wick- 
rotated) time dependence. As we are considering a bosonic 
field that obeys periodic boundary conditions in r, the 
eigenvalues are the known Matsubara frequencies w„i = 
27rm//3Sc, m S Z, and the eigenfunctions are gm{T) = 
exp(— iWrnT). Then, from Eq. ([3]), the coefficients (j)nm 
satisfy the differential equation 

- - [A^ + 0.1] <t^^Us) + Vnmisl (8) 



which can be integrated to give 



4>nm{s) 



(9) 



where the noise coefficients satisfy, from Eq. Q 

(?7nm(s)?7*/„/(s')) = 2kBT5{s - s')5nn'5mm'- (10) 

The formalism developed in fll' allows the average stress 
tensor to be calculated by substituting expression (|6| into 
the stress tensor and taking the average over the fluctua- 
tions, r;(T, r; s), in the limit s — > oo: 



(T(r)) = lim(r[0(s),0(s),r]) 



(11) 



lim y 

S— >C30 ^- — ^ 
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The two functions (j)nm arc replaced by their values in 
terms of the fluctuations, given by Eq. resulting into 
a double integral over two pseudo-times. Then, the aver- 
age (•) is carried out using Eq. (10), that eliminates one 
integral, with the final result: 
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T,m(r) 



/3 ^ A2 + , 



(12) 



where T'„„i(r) = ■!'[/„, /„,!■]■ Moreover, the sum over the 
Matsubara frequencies, can be carried out to obtain 
the resuh 
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1 



(13) 



which can be related with the quantum fluctuation- 
dissipation theorem for the EM field (6|[T6]. Finally, to 
obtain the Casimir force over a certain body, the stress 
tensor must be integrated over the surface fl that defines 
the object 



(T(r)) . dS. 



(14) 



Equation ( 13 ) is the main result of this letter. It gives 



an expression for the quantum Casimir force including the 
effects of a finite, nonvanishing temperature, in terms of 
the eigenvalues and eigenfunctions of the Laplacian op- 
erator. This expression has been easily obtained using 
the stochastic quantization approach to a quantum field, 
together with the formalism of the Langevin equation cal- 
culation of fluctuation- induced forces [TT] . 

Let us note that, as it occurs in most of the calculations 



of the Casimir Force, the expression for (T(r)) in Eq. (13) 



is generally divergent at every point of space r if the sum 
over eigenvalues runs to infinity. However, the expression 
for the force, that is obtained integrating over the surface 
of the body, Eq.(14), is finite. It means that the inte- 



gration over the body regularizes the divergences of the 
averaged stress tensor. If we assume that such regular- 
ization is carried out mode by mode, we can interchange 
the integration over the surface and the summation over 
eigenvalues, to obtain. 



he 
T 



A,i 



1 



i"r„„(r)-rfS. (15) 
Jn 



which is a finite result. Therefore, the interchange of 
the integral and summation regularizes the Casimir force, 
avoiding the use of ultraviolet cutoffs. Other regulariza- 
tions, that in some cases may lead to non-universal forces 
or fluctuations are, for instance, the subtraction of the 
vacuum stress tensor 17 or by averaging the stress tensor 
over a finite area or a finite time [l8|. Having regularized 



the divergences, Eq. ( 15 ), provides a new method to calcu- 



late Casimir forces for a given geometry by diagonalizing 
the Laplace operator. So, this approach is suitable for 
numerical calculations of Casimir forces in complicated, 
realistic geometries. Moreover, this method provides the 
force directly, not as a difference of the free energy with 
respect to a reference state, which in some configurations 
it may be difficult to establish. Also, it can be used as 
a starting point for a perturbative theory for, e.g., non- 
flat geometries, rough surfaces or similar problems. Other 



authors have obtained expressions for the free energy in 
terms of the eigenvalues of the spatial operator, but not 
for the force 
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The expression for the Casimir force, Eq. (13), allows 



one to evaluate the quantum limit, that is, by setting the 
temperature equal to zero (or /3 — > cx)), where the second 



summand inside the brackets in Eq. (13) vanishes, i.e.. 



lhn(T(r))=.^^^""(^) 



(16) 



2 ^ A„ 

In the opposite, classical limit, when /i — > 0, a Taylor 



expansion of the square bracket in Eq. ( 13 ) gives 



(17) 



This expression for the Casimir force has been used 
for classical systems 11 , such as liquid crystals [7j or 



reaction-diffusion systems 21 . The two limits, quantum 



(16 1 and thermal (17), show that the driving force of the 



fluctuations has different origin. In the first case, the pres- 
ence of the factor h indicates the quantum nature of the 
fluctuations, whereas in the second case, the factor /3 re- 
veals its thermal origin. 

Fluctuations. The Casimir force has its origin in fluc- 
tuations, so it is a fluctuating quantity itself. The net 
force is calculated as the average of the fluctuating one, 



Eq.(14). However, the formalism developed in this letter 



allows the calculation of fluctuations of the force, defined 
as customary for the fluctuations: 



/ i'([T(ri).dSi][T(r2)-dS2])-F^ 
Jn Jn 



(18) 



Calculation of ap requires the evaluation of the product 
of 4 fields (/>, as the tensor T is bilinear form of <j). Sub- 
stitution of (f) in the tensor leads to the correlation of a 
product of four noises, ij. Because of the Gaussian nature 
of 77, it factorizes into three products of pairs of noises, as 
shown in Eq. ([5|. The precise expression for the product 
of the two stress tensors in Eq ( 18 ) reads: 



(T(ri)T(r2)) - 



(he) 



4 ^ 



P{Xn)P{Xm)x 



i(ri)r™,„(r2) + 2T„™(ri)T,„„(r2)] , (19) 



where 



P(A) 



(20) 



The first term in ( 19 1 exactly gives when substitued in 



Eq.(18), while the other term (where the factor 2 comes 
because the two last terms in Eq. ^ give equal terms) give 
nonvanishing contributions to a p. Again, there is no need 
to introduce cutoff in the eigenvalues, as interchanging 
the summations with the integral over the bodies regular- 
ize the divergences. Therefore, we find a universal form 
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for the fluctuations, as opposite to other authors [8 18 
The difference has its origin in that we do not compute 
the stress fluctuations or the fluctuations of the force on 
each side of a plate. Rather, we first compute the total 
fiuctuating force on the body (which is finite) and then we 
compute its variance. 




L 



oo 



Fig. 1: Geometry of the considered problem. The Casimir force 
is evaluated for the plate, of area A and general geometry. The 
net force is obtained as the force exerted by a plate at distance 
L and another plate at distance Loo — > oo. The cylinder is 
oriented along the x axis. 

Application to the calculation of the Casimir force of 
a cylinder of arbitrary cross section. In order to show 
how the formalism developed in this letter works, let us 
consider a simple case of a piston depicted in Fig. [T] It 
consists of a piston of area A, but general shape, made of 
a perfectly conducting metal surface 



22 . Two flat con- 



ducting plates of the same cross section of the piston are 
placed at a distance L apart along the x direction. The 
plates are perpendicular to the surface of the cylinder. We 
calculate the Casimir force, Fc, for this plate by evaluat- 



ing the expressions Eq. (131 and later integrating over the 



surface. In order to obtain a finite result, we need a third 
auxiliary plate located at infinite distance, Loo — ^ oo. 

First, we solve the eigenvalue problem for the EM field 
and apply Eqs. (12) and (14) to obtain the force. For the 



EM field the normal component of the stress tensor reads 
TTajK = — ^E^— ^B^, which has to be averaged over 

the noise and then integrated over the surface of the plates, 



as shown in Eq. (14). For perfectly conducting plates, the 
boundary conditions are: E x n = and B • n = 0, where 
n is the normal vector at the surface. In this geometry, the 
EM field can be decomposed into transverse electric (TE) 
and transverse magnetic (TM) modes, which are discussed 
independently [23] . 

For the TM modes, the magnetic field is transversal to 
the X direction, and the vector potential A for TM modes 
can be written as 



Here the fields C{x) and D{y±^ satisfy 



(21) 



vector potential is A = {0, -S d:,N, S dyN)e'''^\ where 



dlC{x) 
Vii?„(r^) 



— k'^C{x), (Neumann BC on x - 
-\lDn{r±), (Dirichlet BC on S), 



where S is the surface of the cylinder and rj^ = iUi^)- 
For the TE set, the electric field is transversal to a;, so the 



the functions S{x) and N{r±) satisfy Eqs. (22) with the 
opposite boundary conditions: Dirichlet for S{x) and Neu- 
mann for iV(rj^). However, in this case, we must exclude 
the constant eigenfunction, with eigenvalue equal to 
zero, as it gives that A = and then E = B = 0. 

Substitution of the TE modes into the expression for 
the stress tensor, and integration over one side of the 
plates gives, after a long but straightforward calculation, 
the Casimir force as 



1 side 



PL 



EE 



(23) 



where fc^ = (na;7r/L)^, and a;„j is defined after Eq. ([t]). 
For the TM modes, one obtains exactly the same expres- 
sion, but are the eigenvalues of the two-dimensional 
(2D) Laplacian with Neumann boundary conditions. We 
will denote the complete set of eigenvalues of the Laplacian 
with Neumann (excluding the zero eigenvalue) and Dirich- 
let boundary conditions by the index p. The expression 



above is the equivalent of Eq. ( 12 ) when the spectrum can 



be split into a longitudinal and transversal part, that is, 
= '^^x + '^p- These series are divergent, but the net 
Casimir force, which is the difference between the force 
exerted from the plate at distance L, and the plate at 
Loo — > oo (see Fig. 1), is finite. 
The sum over the variable Ux 



in Eq. (23) can be car- 



ried out with the help of the Chowla-Selberg summa- 
tion formula 24 . This formula extracts the divergent. 



L-independent part of the summation, which cancels when 



the integral in Eq. ( 23 ) is performed for both sides of the 
piston. Equivalently, the L-independent contribution is 
the same for the second plate at distance L or at distance 
Loo, resulting in 



Fc 



^EE 



/3 



1 



(24) 



Here, A = 2tt/ (3hc is the inverse thermal wavelength. This 
expression gives the finite or regularized Casimir force be- 
tween two plates at distance L, valid for any cross section 
and temperature. The precise geometry of the plates en- 
ters into the double set of eigenvalues of the Laplacian 
(with Neumann and Dirichlet boundary conditions) A^. 

We can proceed to evaluate the Casimir force at T = 0, 
that is, the purely EM case without thermal corrections. 
This case is obtained by noting that, when T — [equiv- 
alent to the limit A — )■ in Eq. (24)], the sum over m 



can be rigorously replaced by an integral. Computing the 
integral for finite A and taking the limit A — 0, the result 



0,L) 
(22) 



Fc 



he 

'2^ 



^^A2jifo(2nLAp)+i^2(2nLAj,)]. (25) 



Here Ka{x) is the modified Bessel function of order a. In 
Ref. 1 19 it was obtained a formula for the free energy of the 
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configuration considered here, that, after differentiation 
with respect to the distance between the plates, leads to 
the force above. 

For short distances, however, the summation above can 
be calculated without explicitly knowing the eigenvalues 
of the Laplacian. Such eigenvalues, Ap, must scale with 
the inverse of the square of the typical size of the pis- 
ton. So, for distances L much smaller than the section of 
the piston, we can replace the sum over the eigenvalues, 
Ap, by an integral, using the asymptotic expression for 
the density of states of the Laplacian in two dimensions: 
p(Ap) = 5^Ap, |25 for each set of Dirichlet or Neumann. 
This result will be universal, independent of the shape of 
the section of the piston, as the density of states is it- 
self universal. The resulting integrals can be performed to 
obtain the Casimir force as 



T=0 



hcTT 

240L4 



(26) 



which is the well-known result of the EM Casimir force for 
infinite parallel plates [l]. In the opposite limit, when L 
is much larger than the typical size of the plate, the argu- 



ment of the Bessel functions in Eq. ( 25 ) is much larger than 
one. Because of the exponential behavior of the Bessel 
functions, only the smallest eigenvalue A^ contributes to 
the sum, with the result 



Ft=q — — - 



he 



2VttL 



51-^1 



3/2 -2LAi 



(27) 



being 171 the degeneration of Ai. Here, a counterintuitive 
result is obtained. One would expect that the thin piston 
would tend to the known one-dimensional (ID) Casimir 
force, but instead an exponential decay of the force is 
found. The known ID result would be obtained if the zero 
eigenvalue were considered. However, this eigenvalue is 
excluded, as it leads to a vanishing eigenfunction. There- 
fore, the ID Casimir force cannot be obtained as the limit 
from three dimensions (3D) to ID. 

At intermediate distances, that is, L comparable to the 
size of the plates, one must solve the eigenvalue problem. 
To illustrate the intermediate behavior and the transition 
from Eq. (26 1 to (27), we study the case of a circular 



cylinder of radius R. In this case, the eigenvalues are the 
zeros of the Bessel function J^(r) and its derivative. The 
dependence with the radius of the cylinder is Xp{R) = 
X^{R — l)/i?^, so the Casimir force as expressed in Eq. 
(25) is a function oi L/R when the force is multiplied by 



R'. 

Figure [2] shows the results for the Casimir force, Eq. 
([25]), when TV = 10,7V = 100 and N = 1,000 eigenval- 
ues are summed (taking into account their degeneration). 
We have also plotted the two limiting results: (i) the 3D 
Casimir force ([26]), valid for L ^ i? with an algebraic 
behavior: FcR^c*: {L/R)~'^ (solid line), (ii) the far dis- 
tance limit, given by (27) with the smallest eigenvalue 
Xi{R) — 3.39 /R^ with degeneracy gi = 2, that is valid for 



L ^ R (dotted line). The transition between both regimes 
is observed at i ~ i?. As expected, when few eigenvalues 
are summed, for instance N = 10, the resulting force is 
only valid in the limit of long distances. As the number of 
eigenvalues increase, the numerical result approaches the 
3D Casimir force. For iV = 1,000 eigenvalues, we have 
excellent results for L/R > 0.05. We remark, however, 
that the full curve can be obtained by only considering 
TV — 100 eigenvalues for large distances and matching this 



numerical result with the asymptotic expression ( 26 ) for 
short distances, with a crossover distance L « 0.3i?. 
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Fig. 2: Casimir force on the plates of a cylindrical piston of 
circular cross section of radius 7? as a function of the distance L 
between the plates. Different curves are obtained by summing 
A'^ eigenvalues of the Laplacian, using Eq. ( 25 I , where A*' is 
indicated in the legend. Solid line is the 3D Casimir force, 
Eq. 1261) and dotted line is the far distance result Eq.(|27l). 



In a similar fashion, we can calculate the thermal 
Casimir force when /i — > 0, or A — > cxd. Then, in Eq. 



(24), only the term with rn = is different from zero, as 

\ ^ A„ 



a result, 



Fh=o — 



p 



1 



(28) 



For short {L<^R) and long (L^R) distances, (28) reads 



F ^^^^ A 



i^fco--^5iAie-2i^^ (29) 



We finish the application to the cylinder by evaluating 
the fiuctuation of the Casimir Force. For the piston geom- 
etry considered in this letter, the fiuctuations of the force 



can be obtained by evaluating Eq. (18). In this case, and 
because of the geometry of the problem, each of the sum- 
mands that appear in the four-point correlations of the 
noise gives F^. Therefore, the fiuctuations of the force, 
for any temperature and cross section, are simply 



2F^. 



(30) 



Somehow similar fiuctuations have been obtained for 
a purely thermal force |8 for each mode that enters in the 
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Casimir force. In a non-equilibrium hydrodynamical sys- 
tem the fluctuations have been measured in (26] by means 
of numerical simulation. The fact that the fluctuation of 
the force is as large as the force itself is a signature of 
fluctuation-induced forces. 

Summary. We have shown in this letter that stochastic 
quantization together with the Langevin formalism pro- 
vides a new approach to calculate Casimir forces in the 
quantum electrodynamics (QED) case, including thermal 
effects. The starting point is the calculation of the Casimir 
force via the stress tensor, which is a function of the fluc- 
tuating fields. Parisi and Wu derived a Langevin equation 
to describe the dynamics of such fields, which can be in- 
tegrated to give an expression for the force. The method 
presented herein is quite simple, and avoids some technical 
complication of other approaches. Moreover, it calculates 
the force directly, instead of the free energy. It also can 
be extended to calculate torques in asymmetric configura- 
tions. Another advantage is that it provides a numerical 
method to calculate the Casimir force in complicated ge- 
ometries, such as those of interest for microelectromechan- 



ical systems (MEMS) devices, alternative to others 27 



The method only requires spectral decomposition of the 
Laplacian operator in the given geometry, and summation 
of the eigenvalues and the integration of the eigenfunction 



along the boundary of the object, as shown in Eqs. (131 
and ( 14 ) . Quantum (T — 0) and classical {h — )■ 0) limits 



are recovered by Eqs. (16) and (17) respectively. The in- 



tegration of the eigenfunction over the surface of the body 
leads to a regularization of the Casimir force, producing 
finite results for the force and the fluctuations. 

This approach can be generalized to include non equilib- 
rium situations. For instance, it can be used to calculate 
the Casimir force between two bodies at different temper- 
Also, the method could be applied to time 



atures 28 29 



evolving temperatures 30 and space dependent tempera- 
ture r(r). 
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